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ABSTRACT 

We present near-infrared (NIR) broadband (0.80-2.42 \xm) spectroscopy of two low 
mass X-ray binaries: V404 Cyg and Cen X-4. One important parameter required in 
the determination of the mass of the compact objects in these systems is the binary 
incUnation. We can determine the inchnation by modehng the elHpsoidal modulations 
of the Roche-lobe filling donor star, but the contamination of the donor star light from 
other components of the binary, particularly the accretion disk, must be taken into 
account. To this end, we determined the donor star contribution to the infrared flux 
by comparing the spectra of V404 Cyg and Cen X-4 to those of various field K-stars of 
known spectral type. For V404 Cyg, we determined that the donor star has a spectral 
type of K3 III. We determined the fractional donor contribution to the NIR flux in the 
H- and K-bands as 0.98 ifc .05 and .97 ± .09, respectively. We remodeled the H-band 



light curve from lSanwal et al.l ( 19961 ) after correcting for the donor star contribution to 



obtain a new value for the binary inclination. From this, we determined the mass of 
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the black hole in V404 Cyg to be Mbh = Q.OlgAfQ. We performed the same spectral 
analysis for Cen X-4 and found the spectral type of the donor star to be in the range K5 
- MIV. The donor star contribution in Cen X-4 is 0.94 it .14 in the H-band while in the 
K-band, the accretion disk can contribute up to 10% of the infrared flux. We remodeled 



the H-band light curve from IShahbaz et al J (11993), 

again correcting for the fractional 



contribution of the donor star to obtain the inclination. From this, we determined the 
mass of the neutron star as M^s = I.SI^Mq. However, the masses obtained for both 
systems should be viewed with some caution since contemporaneous light curve and 
spectral data are required to obtain definitive masses. 

Subject headings: binaries : close - infrared : stars - stars : individual (V404 Cyg, Cen 
X-4) - stars : neutron 



Introduction 



V404 Cyg and Cen X-4 belong to the class of low-mass X-ray binaries (LMXBs) in which 
mass transfer from a late type star to a black hole (BH) or a neutron star (NS) takes place via 
an accretion disk. By determining the mass of the compact object, we can constra in models for 



Frver &: Kalogera 



2001 



Nelemans &: van dan Heuval. 112001 



the formation and evolution of black holes and neutron stars in clo se binary systems (jBrown et al. 
1998 



Of the 16 BH systems with BH 

masses determined to accuracies of 5-30%, most of the masses range from 4-16 Mfn with only one 



syste m, J0422-h32, with a BH mass <5 M© (jcharles k Coell20od : [casares 11200^ : ICehno fc Harrison 



20031 ) . This result is not in agreement with the predictions of BH mass distributions from theoretical 
modeling wh ich predict a continuous decreasing mass d i stribu tion from neutron star masses up to 
10 - 2OM0 (iFryer &: Kalogera I I2OOII : iBelczynski et al.l l2008l ) . However, there currently are not 
enough systems with precise mass determinations to allow us to distinguish between observational 
biases and/or theoretical limitations as the source of the disparity. To increase the number of 
LMXBs with precise compact object mass determinations, we have undertaken NIR spectroscopy 
of several systems to determine the relative contributions of the donor star and the accretion disk 
to the NIR light curves, from which accurate binary inclinations and, in turn, stellar masses are 
derived. In this paper, we present results for two LMXBs, the BH binary V404 Cyg and the NS 
binary Cen X-4. 

V404 Cyg went into outburst on May 22, 1989 and was subsequently discovered with the 
all-sky monitor aboard the Ginga satellite. A high X-ray luminosity of ~ 10^^ ergs/sec suggested 
that the system was a close binary system with a black hole or a neutron star as the accretor 
(|Marsdenlll989l ). During the outburst, the visual magnitude of V404 Cyg increased from V= 18.3 
to V=11.6. One year after its discovery, it had fallen back to V=18.0 ([Wagner et al.lll992l ). Since 
then, V404 Cyg has been studi ed in quiescen c e at o ptical and infrared wavelengths, establishing 
most of the system parameters. ICasares et al.l (119921 ) established the mass function, the minimum 
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mass o f the compact object, as f{M) = 6.26ib0.31MQ which was later refined bv lCasares &l Charles 
(|l994l l to be /(M) = 6.08 ± O.OeM©. They determined the donor star to be of spectral type KOIV 
and obtained the binary mass ratio from the rotational broadening of the absorption lines in the 
donor star spectrum. The ellipsoidal modulation of the donor star has been mode led in the NIR 



to det ermine the bi nary inclinati o n and the black hole mass: in the K-band by IShahbaz et al 



(| 19941 ) and in H by ISanwal et al.l (jl996l ). The resulting black hole mass determinations were 12 
MfT^ and <12.5 Mfn , respectively, where one of the key differences between their results was that 



Shahbaz et al. (ligM l assumed that the donor star was the only sou rce of flux in the NIR , while 
Sanwal et al.l (|l996l ) found evidence for accretion disk contamination. IShahbaz et al.l (|l996l ) tested 
this assumption by performing NIR spectroscopy of V404 Cyg. By comparing the NIR spectrum 
to scaled template spectra from spectral type standard stars, they determined that the donor star 
contributes lOOztll % of the K-band flux. However, for this analysis they relied on fits to the ^^CO 
bands in K, which in other LMXBs (and catac lysmic variables) ha ve been found to be affected by 
metallicity variations in the carbon abundance (jFroning et al.ll2007l ). As a result, a broader spectral 
analysis that examines multiple atomic and molecular species is required to firmly tie down the 
donor star contribution in the NIR. 



Cen X-4 was discovered on May 23, 1969 by the Vela 5B satellite (jConner et al.lll969l ) when 
it went into X-ray outburst. During its second ou tburst in 1979, it incr eas ed in visual magnitude 
by ~6 magnitudes from its pre-outburst phase (jCanizares et al.l 119801 ). iMatsuoka et all (jl980l ) 
identified the mass ac creting star as a neutron star due to the fact that it displayed a Type I burst. 
Shahbaz et al.l (|l993l ) modeled the H-band light curve to obtain a binary incl i nation of 31°-54° 
and the mass of the neutron star as O.5-2.1M0. In quiescence, I Torres et al.l (120021 ) performed 
high-resolution optical spectroscopy on Cen X-4 and determined the mass of the compact object as 
O.49M0 < M < 2.49 Mr7^ . They deterniined the mass function to be /(M) = O.22O±O.OO5M0, which 



was later refined bv lD'Avanzo et al 



(120051 ') to /(M) = 0.201 ± O.OO4M0. The system parameters 



were again refined bv lCasares et al.l (|2007l ) and much tighter constraints were imposed on the radial 
velocity semi-amplitude, K2 = 144.6 ± .3 km s'^, the orbital period, P = .6290522 ± .0000004 d 
and the binary mass ratio, q = .20 ifc .03. The rn ass function adopted for this paper is calculated 
from the P and K2 values of lCasares et al.l ((2003), which is /(M) = 0.197± O.OOIM0. However, as 



with V404 Cyg, in all the above cases, it was assumed that the donor star was the only NIR flux 
source, an assumption that has not been tested in Cen X-4. 

In this paper, we present broadband (0.8-2.42 /xm) NIR spectroscopy of V404 Cyg and Cen 
X-4 obtained with SpeX on the IRTF. By comparing the spectra to field stars of known spectral 
type, we determine the donor star contribution to the infrared flux for both systems. We use these 
results to remodel previously obtained light curve data to determine the binary inclinations and to 
derive precise compact star masses for both targets. 
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2. Observations 

V404 Cyg an d Cen X-4 were ob served on June 7-9, 2007 with SpeX at the NASA Infrared 
Telescope Facihty (jRavner et al. The weather was clear with good seeing conditions (< 0.7"). 

All observations were made using the Short XD (SXD) mode through the 0.5" slit (R=1200), which 
covers the entire wavelength range, 0.8-2. 4/xm, simultaneously in 6 orders. Due to the good seeing 
conditions and the narrow slit used, the observed spectrum is uncontaminated by an accidental 
companion star at 1.5" from V404 Cyg. An AOV star was observed hourly for the purpose of 
telluric correction. Observations were taken in ABBA pairs in two positions along the slit with a 
300 sec exposure time at each position. The slit was oriented along the parallact ic angle throughou t 
the observations, limiting the relative flux calibration uncertainties to < 2% (jVacca et al.ll2003l ). 
Table 1 lists the observations along with the corresponding cumulative exposure times and orbital 
phases observed. 

Data reduction consisted of calibration (sky subtraction, flat-fielding, wavelength calibration), 
optimal spectral extract ion, and telluric cor rection, all of which were carried out using the idl 
based package Spextool (jCushing et al.ll2004l ). Each exposure was shifte d to the rest frame of the 
donor star, u s ing th e orbital ephemeris from ICasares &: Charlea (| 19941 ) for V404 Cyg and from 
Casares et al.l (|2007l ) for Cen X-4. The individual shifted exposures were median-combined and all 
6 orders were merged to obtain the time-averaged spectrum. Using the magnitudes of the AOV star, 
the telluric correction tool, xtellcor, also flux calibrates the spectra. Given our excellent observing 
conditions, the absolute fluxes should be fairly accurate, although we have not made any attempt 
to quantify the accuracy of those values, as they are not necessary for our analysis. 

Figure 1 shows the time-averaged spectrum of V404 Cyg. The spectrum has been boxcar 
smoothed by 3 pixels, corresponding to one resolution element. Statistical errors were propagated 
through all the data reduction and processing steps. The spectrum was corrected for interstella r 
reddening using a dereddening value ofE(B-V)= 1.303 (jHvnes et al.ll2009l : lMiller- Jones et al.ll2009l ). 
The spectrum in red in Figure 1 is the dereddened spectrum of V404 Cyg. Shown in Figure 2, 
3, and 4 are the J-, H-, and K-band spectra. We ha.v e labeled the p r omin e nt spectral features 
using li ne identificat i ons fr om Wallace Sz Hinkle ( 19971 ). IWallace et al.l ( 2000l ) . i Kleinmann Sz Hall 
( 19861 ). iMeyer et al.l ( 19981 ). iKirkpatrick et al.l ( 19931 ). and the Atomic line liso Figure 5 shows 
the time-averaged spectr um of Cen X-4. I t has been boxcar smoothed by 3 pixels and dereddened 
assuming E(B-V)= 0.1 (jBlair et al.l Il984 ). Figures 6, 7, and 8 shows the J-, H-, and K-band 



spectra of Cen X-4. We estimated the S/N in the time-averaged spectra of V404 Cyg and Cen X-4 
by fitting straight lines through various continuum-dominated regions and measuring the scatter 
around these fits. For Cen X-4 we obtained a S/N ~ 15-20 in J and K, and ~ 25-30 in H. As 
the donor in V404 Cyg is likely a giant or sub-giant, there are fewer regions of pure continuum. 
Therefore, we obtained a lower limit to the S/N of ~ 100 in H and K, and ~ 50 in J. 



^http://www.pa.uky.edu/ peter/atomic 
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3. Analysis 
3.1. V404 Cyg 

The spectrum of V404 Cyg shows prominent narrow absorption hnes of neutral metals including 
Al I, Na I, K I, Mg I, Fe I, Si I, Mg I, Ti I and Ca I . The H-and K-bands also show molecular 
absorption bands of ^^CO and ^^CO. The absorption lines originate in the donor star. In contrast 
to many o ther X-ray biiiaries, V404 Cyg does not show any strong emission lines in its infrared 



spectrum. IShahbaz et al.l (|l996l ) performed K-band spectroscopy of V404 Cyg and found no strong 
emission features either, implying that the accretion disk contribution in V404 Cyg in quiescence 
is probably very low at all times. In the K-band, we see both ^^CO and ^'^CO features which are 
a strong indication of an evolved donor star. In this paper, we make use of the entire broadband 
infrared spectrum to ascertain the spectral type of the donor star and then obtain a precise value 
for the relative contributions of the donor star and the accretion disk. 

First, we compared the broadband spectrum of V404 Cyg to field stars of different spectral 
types. These spectra were obtained fro m the IRTF Spectr al Librarj|§ which consists of SpeX SXD 



spectra of stars of known spectral type (jRavner et al.ll2009l ). The template stars used in our analysis 



have roughly solar metallicities and similar S/N to each other. We only used non- variable stars for 
our analysis. Earlier work done on the spect ral type of the donor star in V404 Cyg suggested a 



KOIV spectral type (jCasares fc Charleslll994l ). In the IRTF spectral library, the closest available 
spectral type to a KO IV is a Kl IV. In Figure 9, we compare the spectrum of a Kl IV (IID165438) 
star to the spectrum of V404 Cyg, normalizing the spectrum to that of V404 Cyg just blueward 
of 2.29/xm. It is evident that the CO bands in the template spectrum are much weaker compared 
to those in V404 Cyg. Also, if the template is normalized to V404 Cyg in K, its fiux exceeds the 
J-band flux in V404 Cyg flux by > 30%, indicating that a Kl IV spectrum is too blue to match the 
donor star spectrum in V404 Cyg, unless there is another NIR flux source diluting its contribution. 



Hvnes et al.l (|2009l ) did a multiwavelength study of the spectral energy distribution (SED) of 
V404 Cyg in quiescence and found no substantial contribution from the accretion disk from near- 
UV to near-IR. They also found that for the donor, the spectrum of a KO III star was in agreement 
with the data from UV to near-IR. At the same time, their multiwavelength SED comparison with 
a G8 III and a K2 III star could not be distinguished from a KO III. We find that a K2 III (HD 
132935) spectrum exceeds most of the J-band (.8 - l.l/x m) fiux by >7% when normalized in K, 
indicating that a K2 III donor star spectral type would require a diluting flux source in the NIR. 
Also, the CO bands in a K2 III appear too weak to match that of V404 Cyg in the K-band. In 
Figure 10, we compare a K3 III (HD 221246) spectrum to that of V404 Cyg. In this case, the match 
is much closer, both in the broadband shape and spectral fluxes and in the relative depths of the 
CO bands. When normalized to a 100% in K, we find that the K3 III template exceeds the J-band, 
near the band edge (.8 - .9 m) by at most 5%. On the other hand, the K4 III (HD 207991) 



^http://irtfweb. ifa.hawaii.edu/ spex/WebLibrary /index. html 
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spectrum falls below the V404 Cyg spectrum in J by 20%, so that if that spectral type is correct, 
the donor star may dominate the NIR flux in H and K, with an increasing disk contribution at 
bluer wavelenghs. In conclusion, both the spectral shape and the strength of the CO bands in V404 
Cyg are inconsistent with a Kl IV spectral type for the donor star and indicate that a more evolved 
spectral type is better match to the donor star. In the next section, we calculate the equivalent 
widths of various lines of V404 Cyg and compare that to the same features in field star populations 
to determine the spectral type more precisely. 



3.1.1. Equivalent width analysis 



Forster Schreiberl ([20001), hereafter FS, used the equivalent width (EW) ratios of various NIR 
absorption lines in a population of field giants and supergiants (also a few dwarfs) as indices to 
determine stellar effective temperature and gravity as well as evidence of dilution from non-stellar 
sources. In Table 2, we list the equivalent widths for V404 Cyg for the diagnostic features in FS. The 
error bars represent uncertainties in continuum placement only. By comparing our EWs to those in 
Figure 5 of FS, we obtain a donor effective temperature range of 4200-4600 K, implying a spectral 
type between K2 III and K4 III, assuming that the donor star is the only NIR flux source. FS also 
demonstrated how the EW ratios of certain lines would change if there was a source of diluting 
continuum emission. In their Figure 8 they plotted dilution- free indices (ratios of atomic/molecular 
EWs) and directional indicators to show how the indices would change as the contribution from a 
diluting source increases. A comparison of the same EW ratios in V404 Cyg to their indices shows 
that the V404 Cyg ratios fall on their plots for the non-diluted sources, suggesting that the amount 
of dilution in V404 Cyg in both H and K is negligible. 



3.1.2. Field star fits to V404 Cyg 

To obtain the fraction of NIR continuum emission originating from the donor star, /, we 
compared the spectrum of V404 Cyg with those of field stars ranging from K2 III to K4 III spectral 
types. The spectrum of the template star was convolved with a Gaussian of 38.8 km s~^ FWHM 
in order to take into account the rotational velocity of the donor star in V404 Cyg (but we note 
that this correction was negligible given the 250 km resolution of our spectra). Before fitting, 
we normalized the spectra of both V404 Cyg and the template by fitting cubic splines along the 
local continuum surrounding the lines under analysis. The continuum points were picked by eye. 

We determined the donor fraction as follows: after normalizing both the target and the tem- 
plate spectra, we multiplied the latter by a fraction, /, which was incremented from to 1 in steps 
of 0.01. The result was subtracted from the spectrum of V404 Cyg to obtain the residual. We 
then minimized the chi-square, xt^ between the residual and the mean of the residual to obtain the 
donor fraction that gave us the best fit. In cases where the template had weaker absorption features 
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compared to V404 Cyg, we let the donor fraction, /, take on values greater than 1. Even though 
this represents an unphysical scenario, it is useful as a guide to constrain the correct spectral type 
of the donor star by calculating the extent of weakness of the absorption feature in a particular 
template. For the fitting procedure, we selected only the strongest absorption features and fit 
absorption features for individual species separately. There were cases where absorption features 
of two different species occurred very close together in wavelength. If we were not able to separate 
them and identify the continuum around those absorption features, we did not include them in our 
analysis. 

In order to better account for the uncertainty associated with the data, we divided our indi- 
vidual exposures into two sets. They were then separately combined and the result of one divided 
by the other. The deviation from the mean was then determined for each spectral order separately 
to assign the errors. Since the J-band comprises of three orders, we obtained the mean of these to 
assign the error bars for the J-band. We used these errors for the analysis for V404 Cyg as it was 
clear that the statistical error bars are too small and do not account for systematic uncertainties 
in the data reduction process. For the template stars we used the statistical errors provided with 
the template spectra, but we note that the uncertainties for V404 Cyg dominate over the error in 
the template stars while fitting various spectral features. 

Table 3 shows the donor fractions calculated for K2 III, K3 III and K4 III spectral types along 
with the Xu values. When we compare the xt values for fits across the entire K and H-band at 
once, we obtain K3 III as the best match for the donor in V404 Cyg. However, in order to make 
our final choice on the donor spectral type, we also looked at all the individual absorption features 
selected for our analysis. 

In the K-band, we fit the CO bands together from 2.29 - 2.42 /xm and obtained the best fit 
with a K3 III spectral type. There is a Na I feature near 2.34 /um included in the fit region but 
the fit is dominated by the prominent CO-features. Generally, the K3 III star has the lowest Xu 
for fits to individual absorption features compared to the K2 III or the K4 III template stars. For 
the K2 III star all the absorption features are consistently weak in the template compared to V404 
Cyg (resulting in the unphysical fit, / > 1). The K4 III template has fewer / > 1 deviations than 
for the K2 III star, but some of the lines cannot be fit with / < 1 and the Xu values are generally 
larger than for a K3 III template. Even for the K3 III template, some lines give / > 1 including the 
Mg I, Ca I and the Al I feature near 2.28 /.im, 1.98 fiui and 2.112 - 2.128/im, respectively. These 
features deviate less from / = 1 for a K3 III star than for the other templates but may indicate 
systematic mismatches (in temperature, gravity, or metallicity) between V404 Cyg and the field 
star templates. Overall, a K3 III star provides the most consistent match to the K-band spectrum 
of V404 Cyg. 

In the H-band, looking at the CO features, we find / = 1.24 for the K2 HI template, implying 
very weak ^^CO features in the template compared to V404 Cyg, whereas in K3 HI and K4 III, 
the value of / is 0.98 and 0.95, respectively. In case of the K2 HI star, the value of donor fraction 
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is consistently > 1 for all the individual features in the H-band also indicating that most of its 
absorption features are weaker than in V404 Cyg. With the K3 III and K4 III stars, we find that 
Mg I is always weaker in the template than in V404 Cyg (except for one feature at 1.71//m when 
compared against a K4 III template). This trend is also observed for Mg I line in K. The same is 
true for the Al I feature in the H-band though it is close to / = 1 for the K3 III template. In the 
H-band too, we find that a K3 III star is the most consistent match to the donor star spectrum in 
V404 Cyg in terms of Xu well as the strength of the lines indicated by the donor fraction. 

In the J-band too, we find that when we compare the Xu values and also the strength of the 
absorption features of the template with respect to V404 Cyg for the three different spectral types, 
a K3 HI spectral type best matches the donor type in V404 Cyg. From our analysis with the 
spectral energy distribution, equivalent width calculation of various lines and direct determination 
of donor fractions, we conclude that the donor star in V404 Cyg has a K3 HI spectral type. We 
see from Table 3 that the Mg I features for a K3 HI spectral type has been consistently weaker in 
the template than in V404 Cyg (except the one feature near .856/im). This suggests an abundance 
mismatch between the template and V404 Cyg for this species. The same is perhaps true for Al I 
but there are not enough features to establish a consistent trend. For this reason, we did not include 
Mg I in our final donor fraction calculation. 

Even though there is a good qualitative match between the spectrum and the K3 HI template, 
our x^ statistics are poor, ranging from << 1 for fits to individual lines in J and H to >> 1 
for fits to individual lines in H and broadband fits in H and K. We attempted to take into account 
the systematic uncertainties that dominate in NIR ground-based spectra by recalculating the error 
bars as described above, but it is clear that our error bars are still too large in J and portions of K 
and too small in H. Due to this, we did not rely on x^ statistics to obtain our final donor fraction. 
Instead, we determined the donor fraction in V404 Cyg by taking an average of the best fits to 
multiple absorption features using a K3 III template and used the standard deviation between these 
line fits to represent the uncertainty on the final value of /. We believe this is a better tracer of 
systematic uncertainties in our fitting process particularly because it also takes into account the 
possibility of slight mismatches between the V404 Cyg and the stellar template. Using this method, 
we obtained a H-band donor fraction oi f = .98 ± .05 and a K-band fraction oi f = .97 ± .09 . 
Thus, we find no subst antial NIR contribut ion from the accretion disk in V404 Cyg similar to 



Hvnes et al.l (12009|) and lShahbaz et al.l (Il99fil l 



In Figure 11, we show the the K-band spectrum of V404 Cyg overplotted with the spectrum 
of a K3 III star scaled by a factor of 0.97. The close qualitative match between the template and 
spectrum is evident. The Mg I lines at 2.10 and 2.28 //m are weaker in the template spectrum than 
in the V404 Cyg. Shown in Figure 12 is the spectrum of V404 Cyg in the H-band over plotted with 
a scaled K3 III spectrum using / = 0.98. The Mg I lines at 1.50^m, 1.58/im and 1.71/um in the 
template are weaker than in V404 Cyg but the other absorption features give satisfactory fits. 
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3.2. Cen X-4 

The spectrum of Cen X-4 shows broad emission hnes of H I and He I and narrow absorption hnes 
of neutral metals including transitions of Na I, Mg I, Al I, Si I, Ca I, Fe I and molecules including 
^^CO. The abs orption lines orig i nate i n the cool donor star while the emission lines originate in the 



accretion disk. IChevalier et alJ (|l989l ) observed Cen X-4 in quiescence and determined the donor 
star to be of spectral type K5V - K7V. They also estimated that 25-30% of the total light was 
contributed by the accretion disk at visual wavelengths. In the infrared, t he contribution of the 



accretion disk to the total flux has traditionally been expected to be even less. IShahbaz et al.l (|l993l ) 
obtained the H-band light curve of Cen X-4 and predicted the contribution from the accretion disk 
to be significantly less than 10%. 

Our spectrum of Cen X-4 is noisier compared to that of V404 Cyg which will make a precise 
spectral type determination more challenging. We compared the spectrum of Cen X-4 with field 
dwarfs ranging from K5 - Ml V. The normalization was done at a wavelength just blueward of the 
^^CO bandhead near 2.29 //m. If the donor star in Cen X-4 is the only contributor to the infrared 
flux, we find that the spectral types K5 - K7 V can be excluded. When the flux of a K5V star is 
normalized to that of Cen X-4 in the K-band, it exceeds the flux in the H-band by > 15%. This 
excess in flux decreases as we move to cooler stars. For example, in Figure 13 we find an excess 
of ~ 5% in a K7 V star in the H-band and in the short J-band when normalized in K. From the 
comparisons, we conclude that if the donor star in Cen X-4 is the only NIR flux contributor, then 
it must be of spectral type later than K7 V; otherwise, some dilution must be present. For later 
spectral types, the flux of the donor star falls below that of Cen X-4 in the J-band and matches 
the H-band, as seen for an MO V star in Figure 14. We also find ~ 1% (less than the error in the 
K-band) excess in the shorter K-band in all spectral types including an MO V, indicating that some 
dilution, if quite small, may be needed for later spectral types as well. 



3.2.1. Equivalent width analysis 

The analysis in FS focuses more on giants and supergiants than dwarfs, making use of his 
temperature and dilution indices ideal for V404 Cyg but less so for Cen X-4. However, we calculated 
the equivalent widths of the absorption features used by FS. Table 4 lists the values obtained for 
the equivalent widths of several absorption lines in Cen X-4. Due to the lack of data on dwarfs in 
FS coupled with a low S/N Cen X-4 spectrum, we obtained a large scatter in the T^ff values. We 
relied on just the prominent lines in Cen X-4 like the Ca I line at 2.26/im, the ^^CO bandhead at 
2.29/im, the Na I line at 2.2fim., and the ^^CO bandhead at 1.62/im to constrain the spectral range 
of the donor. For Ca I at 2.26/im in FS, there are no data on dwarfs for the EW we obtained. With 
the remaining lines, we find K7 V - Ml V as the spectral range of the donor star in Cen X-4 in 
the event of no dilution. 
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Ali et al.l (|l995l ) have performed a similar analysis as FS but with emphasis on dwarfs ranging 
from F3 - M6 V. We calculated the equivalent widths for Ca I, Mg I, Na I and the ^^CO (2-0) 
band using the integr a tion l imits defined in their paper, which we show in Table 5. (Note that 
the FS and I Ali et alj (|l995l ) EWs are typically different for the same spectral features because 
of different choices in the integration limits over which t he EWs were ca, lculated.) We compared 
the EWs obtained for Cen X-4 with Figures 3 and 4 of I Ali et alJ (jl995l ). Again, if we use only 
the prominent lines we find that the EW from Na I points to a spectral type K5 - K7 V and that 
Ca I predicts a spectral type later than K7 V. The EW of the ^^CO bandhead at 2.29//m give s 
inconclusive results as it falls within the flat portion of their data as well as Mg 1. 1 Ali et alJ (|l995l ) 
also noted that the sensitivity of CO to gravity makes it less useful as a temperature diagnostic. 
In the absence of dilution, we estimate the range of the spectral type of the donor in Cen X-4 as 
K5 - Ml V. 



Our conclusion from the EW calculation using the absorption features of I Ali et al.l (|l995l ) is 
different from the estimate of FS in that we cannot rule out K5 V as a possible donor type. At the 
same time, from a spectral energy distribution comparison, if a spectral type as early as K5 V or 
K7 V were to be a possible donor in Cen X-4, we should expect to find some amount of dilution. 
We will examine the full K5 V - Ml V range to find out the best spectral type for the donor in 
Cen X-4. In the next section, we present the results obtained from performing a dilution analysis.. 



3.2.2. Field star fits to Cen X-4 

We compared the spectrum of Cen X-4 with those of field dwarfs: K5 V (HD36003), K7 V 
(HD237903) and Ml V (HD19305) spectral types from the IRTF spectral library using the same 
procedure used in calculating the donor star fraction in V404 Cyg. The spectrum of the template 
star was convolved with a Gaussian of 43 km s~^ FWHM to account for the rotational velocity of 
the donor in Cen X-4. The same procedure for determining the donor fraction in V404 Cyg was 
used for Cen X-4. For the fitting process, we chose only the most prominent absorption features 
in the spectrum of Cen X-4. The absorption features of individual species were fit separately and 
when such features occurred too close in wavelength such that the continuum region around an 
individual species could not be separated, we did not include that in our analysis. Due to the lower 
S/N, we fitted fewer lines than in V404 Cyg. 

In order to better account for the systematic uncertainties, we determined the error bars in 
Cen X-4 by fitting straight lines through various continuum-dominated regions in the J-, H- and 
K-bands. For Cen X-4, we find this a more reliable method for error determination than the 
method of dividing the data into half used for V404 Cyg. (Since the NIR spectra of giants have 
fewer regions of continuum without strong line absorption, this method could not be applied in the 
error determination of V404 Cyg.) We propagated the error from the template stars as well in our 
analysis although the xt values are dominated by the uncertainty from the Cen X-4 spectrum. 
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Table 6 shows the best fit donor fractions along with the xt values obtained in Cen X-4 for 
various wavelength ranges with K5 V - MO V template stars. When we compare the Xu values for 
fits across the entire K and H-band at once (after masking out the emission lines), we obtain an 
MO V and a K5 V as the best matches, respectively. The entire J-band could not be normalized 
satisfactorily at once and hence is absent in Table 6. We also looked at the Xu values and the donor 
fractions for individual lines. In the K-band, a K5 V star has the lowest Xu for most of the features 
while in the H- and the J-bands, it is an MO V and a K5 V star, respectively. For all the features, 
we also examined the donor fractions for cases where / > 1. From Table 6, it is seen that in the 
K-band, an MO V satisfies / < 1 for most cases and a K5 V does not; whereas in the H-band, all 
three spectral types are almost equally good, though the Si I feature (1.585 - 1.605 /xm) gives / > 1 
for all three templates. The same can be said of the J-band considering that the fractions / = 1.03 
and / = 1.06 obtained for a K7 V and an MO V in the Mg I feature are close to / = 1 (especially 
given the low S/N of Cen X-4). From this examination, it is evident that that data will not support 
narrowing the spectral type of the donor star from the K5V - MOV range under consideration. 

Three features (the Na I line in the K-band, the Si I features in the H- and J-band and the 
Fe I feature in the J-band) are consistently weaker in all the templates than in Cen X-4. This 
might represent abundance mismatch between the donor in Cen X-4 and the templates akin to the 
mismatch seen in the Mg I feature in V404 Cyg. FS has noted that the temperature and luminosity 
variations of Ca I, Na I and Fe I features in the K-band are primarily governed by a blend of those 
species with Sc I, Ti I and V I which could explain why it is challenging to self-consistently match 
templates to these features. (A similar problem was encountered in K-band fits to the spectrum of 
A0620-00; Froning et al. 2007.) These features cannot be eliminated from the final donor fraction 
calculation, however. Owing to the presence of only a couple of these species in the entire spectrum, 
it is not possible to make a decisive statement about abundance mismatch for specific elements. 
Other notable features are the ^^CO features in the K-band and in the H-band, which in other 
X-ray binaries such as A0620-00 have been anomalously weak. The CO feature in the K-band is 
much weaker in the template for a K5 V and a K7 V compared to Cen X-4, but all of the spectral 
types give consistent results of / ~ 0.93 for the H-band feature. A better S/N spectrum of Cen X-4 
is required to make a conclusive argument about this discrepancy especially as the S /N deteriorates 
sharply in the long K region where the CO bandheads occur. 

We determined the final donor fraction by averaging over all the individual selected features 
from Table 6 for all three of the spectral type templates from K5 ~ MO V. We obtain a mean donor 
fraction of / = 0.94 it .14 in the H-band and obtain / = 1.09 it .20 in K. Thus, the accretion disk 
can contribute up to 10% in the K-band. If a K5 V were to be excluded from the donor spectral 
type (based on its weak K-band features), we find a disk contribution of up to 17% in K. This is 
consistent with the H -band fraction of 0.94, which would imply a ~6 - 10% dilution in K (depending 
on the template spectral type). The H-band donor fraction contribution of / = 0.94 ± .14 remains 
the same whether the K5V star is included o r not. Since the mass of the NS in Cen X-4 will be 
obtained from the H-band light curve data of IShahbaz et al.l (|l993l ). the inclusion or exclusion of 
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a K5 V star will not affect the final mass. Figure 15 and 16 shows a comparison of the spectrum 
of Cen X-4 in the H- and K-band respectively, after scaling a K7 V and an MO V template by a 
donor fraction of .94 and 1.0 respectively. 



4. Discussion 



4.1. V404 Cyg 



4-1.1. Accretion disk contamination and mass of the black hole in V4O4 Cyg 



Previous estimates have given KO IV as the donor type in V404 Cyg (jCasares &: Charleslll994l ). 
Our high S/N NIR spectrum shows strong ^^CO and ^'^CO features that point to a more evolved 
star. By comparing the spectral energy distribution with various field stars, EW diagnostics and 
dilution analysis, we found that a K3 III star is the best match to the donor star of V404 Cyg with 
minimal accretion disk contribution. From a comparison of the spectrum of V404 Cyg with a scaled 
K3 III star, one can appreciate the close qualitative match in the two spectra. In the K-band, the 
CO bands and most metals are consistent with the K3 III spectral type but M g I appears supersola r 



in both the H- a nd K-bands. Unlike other compact binaries like A06 20-00 (IFroning et al 



XTE J1118+480 (IHaswell et al.ll2002l ) and some cataclysmic variables (IHarrison et al.ll200d . l2004l ) 



200711 



that show weak CO bandheads and depleted carbon abundances, the CO bandheads in V404 Cyg 
are consistent with the solar abundance K3III template spectrum. 



Sanwal et al.l (|l996l ) obtained an H-band light curve of V404 Cyg. We retrieved their data from 
their Table 2 and remodeled the light curve taking into account the (small) dilu tion of the donor star 
flux. The light curve modeling program is the same as used for A0620-00 (jFroning &: Robinson 
2OOII ). The parameters used for the model are Tpy f =4300 K, a gravity darke ning exponent of 
0.08 (Lucy, 1967) and a limb darkening coefficient fromlClaret &: Gimenea (119951) fo r a log 5 value 
suitable for a K3 III star (IBerdyugina Sz Savanovlll994l : ICayrel De Strobel et al.lll980l ) . For the error 



bars o n the light curve, we used the scatter on the reference field star observed by I Sanwal et al 
(|l996l ) (.02 mag) and equally weigh ted the points. For a mass ratio of q = 0.060j;QQ5 and /(M) = 



6.08 ± .06 (jCasares &: Charleslll994l ). the mass of the BH can be written to one unknown as Mbh 



= (6.83 zb .09) sin~^ i, where i is the inclination of the binary. We varied the donor fraction from 
.93 - 1.00 and obtained a best fit (x^ = 3.3) to the H-band light curve data for an inclination of 
i = {67~^f)° where the uncertainty in the inclination is propagated from the fractional donor star 
contribution. Figure 17 shows the best fit light curve of V404 Cyg obtained for q = .06 and i = 67°. 
The uncertainty was dominated by the change in the donor fraction. Accordingly, we obtained the 
mass of the black hole in V404 Cyg as 9.O+;gM0 for q = 0.06. 



It has been shown in ICantrell et al.l (|2008l . l20ld ) that fitting ellipsoidal models to the light 
curves of A0620-00 in its passive state increases the binary inclination considerably and results 
in reduced mass of the black hole than previous studies. The variability in the light curve shape 
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(which can change from night to night) has been attributed to the disk component. Since V404 Cyg 
has a high S/N, we tested the stabihty of the donor fraction from the spectrum obtained for each 
individual night and found it to be fairly cons tant for each night . V404 Cyg exhibits fast strong 
variability which is likely produced by the disk (jSanwal et al.lll996l ). To test what effect that would 
have on the inclination and the black hole mass, we picked points from the lower envelope of their 
light curve data in Figure 1 by eye. We found a best fit for an inclination of ~ 72° wh ich translates 
to a r eduction of black hole mass by about IMq which is still in agreement with ISanwal et al 



(1199a). 



Shahbaz et al.l (|l996l ) obtained a 100% contribution of the donor star to the infrared flux 



in the K-band which implies a binary inclination of 56°. Despite our comparable result for the 
minimal accretion disk contribution, their i nclination does not ag ree with our value within the 
error bars. They u sed the light curve from IShahbaz et al.l (119941 ) to determine the inclination. 
Sanwal et al.l (|l996l ) recalculated their K-band model light curve by equally weighing each night of 
data and obtained asignificant improvement on the variance of the model light curve. From Figure 
2 of ISanwal et al.l (jl996l ). for a mass ratio of q=.07, the absolute lower limit to the inclination 
corresponds to ~ 60°. They established an upper limit of 12. SM© for the mass of the black hole 
in V404 Cyg and our result is consistent with their limits with respect to both the inclination and 
the BH mass. 



4- 1-2. Accretion disk contamination and Mass of the neutron star in Cen X-4 



From a comparison of the spectrum of Cen X-4 to those of various field dwarfs, we found that 
the donor star in Cen X-4 is probably not the only source of NIR flux. We found that in the H-band, 
a fraction, 0.94 it 0.14, of the NIR flux comes from the donor star while in the K-ba nd this fraction 



is 1.0 9 ± .20. Earlier estim ates for the donor type in Cen X-4 include a K7 V star (IShahbaz et al 



19931) , a K3 V-K5V star (ITorres et al.ll2002l l and a K3 V-K5 V star (jD'Avanzo et al.ll200,5l V We 
have estimated the donor type in Cen X-4 to lie between a K5 V and an Ml V. If the donor star 
is a K5 or K7V, there must be at least 5 - 15% K-band dilution to reconcile the NIR spectral 
energy distribution in Cen X-4 with the template SED. A later type donor star of MOV or MIV 
requires only a few percent dilution in K to match the broadband SED. Higher signal to noise data 
is required to constrain the donor spectral type even further. 



We use the orbital parameters from ICasares et al.l (120071 ) and determine the mass function 
is found to be /(M) = .197 it .001. Using this value of the mass function and the binary mass 
ratio, q = .20 ± .03, the mass of the neutron star in Cen X-4 can be written upto one uknown as 
(0.28±.01) sin-^i Mq. 

For modeling purposes, we chose Tg// = 4300 K for a K7 V. We chose a val ue of log q=4.5 ap- 
prop riate for a K7 V star to obtain the corresponding limb darkening coefficients (jClaret Sz Gimenez 
19951 ). We also chose a value of 0.08 for the gravity darkening coefficient assuming that the donor 
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has a convective envelope (jLucvlll967l : ISarnalll989l ) 



The flux from the accretion disk wih reduce the apparent amphtude of the ehipsoidal modu- 
lations which in turn will underestimate the binary inclination and overestimate the mass of the 
neutron star. We determined that the contamination from the accretion disk could be up to 20% 
in H. We can make an initial estimate of the mass of the neutron star in Cen X-4 by combining 
our dilution estimate with the H-band light curve observations of Shahbaz et al. We read in their 
light curve data from their Figure 1. Their data points include uncertainties in both magnitude 
and orbital phase, which was derived from how they binned their data (rescaling their errors in 
each phase bin to get = 1 and then averaging the points). We could not read these values off 
their graph, so we assumed an equal weighting on each point using the 0.006 mag error for their 
reference star measurement. 



In order to take into account the uncertainties in the orbital phase bv lShahbaz et al.l (|l993l ). we 
varied the donor fraction in Cen X-4 from 0.80 < / < 1.00 along with a phase-offset that gave the 
lowest for a certain inclination. We obtained an inclination of (35^^)° for the binary from which 
we obtained the mass of the neutron star as M^r^ = 1.5^4-^0. Figure 18 shows the best fit (x^ = 
9.57) to the H-band light curve (dashed line) of Cen X-4 for q = 0.20 and i= 36° for a phase-offset 
of .03, after it was modeled with a K7 V donor plus 6% constant dilution contribution. The solid 
line represents the best fit for no phase-offset (x^ = 12.2). Note, however, that the uncertainty on 
our final result only refiects the uncertainty in varying the donor star fraction. It does not include 
the uncertainty in the light curve measurements. As a result, our derived neutron star mass and 
uncertainty should be viewed with some caution. 

With that caveat, t h e valu e I.IM^ < Mjy . s < I. 6 M0 we find is consistent with past estimates 



found by ICasares et al.l (|2007l ). [Torres et al.l (|2002l ). lO'Avanzo et all (|2005l ) and IShahbaz et al 
(jl993l ). To include the spectral types between MO V - K5 V we varied the temperature, Teff 
from 3700 K to 4500 K, and found that the inclination changes by ~ 1°. We obtain the same 
variation in the inclination if we use logg(=4.0 if the donor in Cen X-4 were to be a sub-giant. The 
mass ratio value also has a negligible effect on the inclination. Upon varying the mass ratio within 
its uncertainties (0.20±0.03), we found the inclination to change by less than a degree. The biggest 
uncertainty comes from varying the donor fraction compared to any other parameters. 



Both Chevalier et al. (jl989l ) and Shahbaz et al. ( 1993 ) have argued that the donor in Cen X-4 
is more likely a subgiant than a dwarf. They pointed out that based on the orbital period alone, 
if the donor were to be a main sequence star, it would underfill the Roche-lobe by a factor of ~ 2. 
The CO bands would provide the best tracer of whether the donor star is evolved. We could not 
test this hypothesis here, both because of the relatively poor data quality in the long wavelength 
end of the K-band where the strongest CO bandheads occur and because of the dearth of subgiant 
stars in the IRTF spectral library. Hence, the possibility that the donor in Cen X-4 could be a 
subgiant (later than Kl IV) remains open at this point . 



It has been seen that pulsar masses are typically higher (~ 1.50 Mq) in neutron star-white 
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dwarf systems th an in neutron sta r-neutron star systems and thi s is attributed t o high er accretion 
in these systems (jNice et al.ll2008l ). It has been pointed out by IShahbaz et al.l (jl993l ) that either 
the neutron star (NS) in Cen X-4 was formed as it reached the end of nuclear evolution in the 
progenitor's core or it could have collapsed to form a neutron star by gradual accretion on to a 
dwarf with electron degenerate core. In the latter case, mass of the neutron star could be higher 
than typical n eutron star masses as they can accrete up to severa l tenths of solar mass during 
their lifetimes (jCasares II2006I : iBhattacharya &: van dan Heuve]|[l99ll ). Our value of the mass of the 
neutron star in Cen X-4 exhibits a range I.IMq < Mj^s < I-SMq with the most probable value of 
1.5Mq. The error on the mass needs to be constrained further in order to realize whether Cen X-4 
could be a potential system for understanding nuclear matter and ruling out diffe rent equations of 
state similar to the neutron star LMXB candidate, Cyg X-2 (jCasares et al.lll998l ). 



In the analysis done for both V404 Cyg and Cen X-4, there are areas that contribute to 
systematic uncertainties which may not be accounted for in our final uncertainties. One of these is 
using a field star to duplicate the Roche-lobe filling donor stars. In V404 Cyg, we found that a K3 
III field star with near-solar abundance matched the donor of V404 Cyg quite well except for the 
Mg I lines. For Cen X-4, we find that after applying the final donor fraction, some of the lines in the 
template were still weaker than Cen X-4. Using a field dwarf with solar abundances to represent 
a Roche-lobe filled star in Cen X-4 could be a plausible explanation for this discrepancy, as could 
minor mismatches between template and donor star in temperature, gravity, or metallicity. We 
attempted to mitigate this potential error source in Cen X-4 by averaging over fits to multiple lines 
and multiple spectral types to obtain the fraction and its uncertainty. Fortunately, our results are 
dominated by the H-band donor fraction where we obtained consistent results for all three template 
stars. 

Previous studies have assumed that the light curve shape is consistent throughout quiescence. 
However, it has been shown in several systems that this assumption is not valid and w i ll lea d 
to systematic errors in the determination of the compact star masses (jCantrell et al.l |2008| . 12010 ) . 
We obtained the spectra of our targets in 2007 whereas the light curve data for V404 Cyg and 
Cen X-4 were obtained in 1993 and 1990-91 respectively. This wide time gap between obtaining 
the light curve data and the spectra adds to the uncertainty associated with the mass of the 
compact objects. In V404 Cyg, the accretion disk contribution is v ery small and has remained so 

Our fractional contribution 



over multiple observations (IShahbaz et al.l 119961 : iHvnes et al 



200 



for the donor star agrees with the value obtained bv lShahbaz et al.l ( 19961 ) within their statistical 
uncertainties implying that the small disk contamination in the K-band has changed negligibly in 
the time gap between the two observations. For Cen X-4, we do not have as much information 
about the long-term variability of the system. However, the uncertainty on our donor star fraction 
(±14% in H) and resulting uncertainty on the neutron star mass is broad enough to span the 
observed range of variability in the system. We also note that the 1.5Mq result is consistent with 
the masses measured in other NS LMXBs. Contemporaneous light curve and spectral data are 
recommended to obtain a definitive NS mass in Cen X-4. 
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5. Conclusion 

We have obtained broadband near infrared spectroscopy of V404 Cyg and Cen X-4 in order 
to directly evaluate the relative contributions in the infrared light due to the donor star and the 
accretion disk and determine precise compact star masses. We determined the spectral type of 
the donor star in V404 Cyg to be a K3 111. We performed dilution analysis based on individual 
prominent metal lines and groups of lines and determined the donor star fraction to be 0.98 =b .05 
in the H-band and 0.97 it .09 in the K-band. With previous light curve data in the H-band, we 
obtained the inclination of the binary as i = (67^'^)° and a mass for the BH of 9.0^:6 Mq . In Cen 
X-4 we performed a similar analysis but restricted our fits to the strongest absorption lines due to 
the lower S/N of the spectrum. We determined the spectral type of the donor star to lie within 
the range K5 - Ml V but could not constrain it further. The dilution fraction was determined to 
be 0.94 ± .14 in the H-band and 1.09 it .20 in the K-band. From past H-band data we modeled the 
light curve with a donor star plus a constant flux from the accretion disk and determined the mass 
of the neutron star in Cen X-4 as 1.5^'4M0 which agrees well with past estimates. These mass 
measurements should be viewed with some caution however, since contemporaneous light curve and 
spectral data are required to obtain definitive masses. 

We would like to thank Mattew Beasley for his assistance in observing V404 Cyg and Cen X-4 
and the IRTF staff for their help. We would also like to thank Niall Gaffney for useful discussions. 
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Table 1. Observations 



Object 


Date 


Exposure time(min) 


Orbital Phase 


V404Cyg 


06/07/07 


100 


0.59-0.61 


V404Cyg 


06/08/07 


190 


0.74-0.76 


V404Cyg 


06/09/07 


150 


0.89-0.91 


CenX - 4 


06/07/07 


100 


0.80-0.92 


CenX - 4 


06/08/07 


30 


a 


CenX - 4 


06/09/07 


110 


0.97-0.05 



''Data from 06/08 for Cen X-4 was not used due to poor data 
quality at high extinction. 
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Table 2. Equivalent widths for V404 Cyg using FS indices 



Feature 


Symbol 


Integration limits 


Equivalent 


width(A) 


Sil 1.5892 


Wi.59 


1.5870-1.5910 


3.36 ± 


0.07 


1^00(6,3) 1.6187 


Wl.62 


1.6175-1.6220 


4.10 ± 


0.08 


Nal 2.2076 




2.2053-2.2101 


2.03 ± 


0.10 


Fel 2.2263 


Wpel 


2.2248-2.2293 


1.22 ± 


0.74 


Fel 2.2387 


WFe2 


2.2367-2.2402 


0.76 ± 


0.12 


Cal 2.2636 


Wca 


2.2611-2.2662 


2.06 ± 


0.23 


Mgl 2.2814 


WMg 


2.2788-2.2840 


1.18 ± 


0.18 


1^00(2,0) 2.2935 


W2.29 


2.2924-2.2977 


9.11 ± 


0.76 


i2cO(3,l) 2.3227 


W2.32 


2.3218-2.3272 


9.15 ± 


0.28 


^^CO(2,0) 2.3448 


W2.34 


2.3436-2.3491 


3.30 ± 


0.56 



No te. — The equiva l ent w idths were calculated using the integration limits 
from iForster Schreibej ||2000|) . The error bar on the measurements are based 
on continuum placement uncertainties. 
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Table 3. Donor fractions, /, obtained for V404 Cyg with K2 III - K4 III field stars 



Band 


Wavelength range(/^m) 


Feature 


/(K2 III) 




/(K3 III) 




/(K4 III) 






1.97-2.42 


entire K-band 


1.05 


8.90 


0.95 


5.83 


0.84 


8.79 




1.97- 1.99 


Ca I 


1.36 


2.55 


1.07 


4.50 


1.18 


9.53 




2.112 - 2.128 


Al I 


1.50 


0.20 


1.05 


0.52 


1.57 


.35 


K 


2.203 - 2.214 


Na I 


1.50 


0.24 


0.94 


0.15 


1.07 


0.23 




2.26 - 2.27 


Ca I 


1.24 


0.21 


0.85 


0.15 


0.98 


0.26 




2.276-2.286 


Mgl 


1.60 


0.89 


1.25 


0.87 


1.87 


0.73 




2.28 - 2.42 


CO bands 


1.05 


13.0 


0.95 


6.70 


0.87 


17.80 



H 



1.48 - 


1.72 


entire H-band 


1.33 


21.2 


1.09 


13.42 


1.05 


18.76 


1.48 - 


1.52 


Mgl 


1.53 


3.57 


1.19 


4.34 


1.40 


5.96 


1.574- 


1.587 


Mgl 


1.43 


7.00 


1.15 


4.99 


1.19 


6.78 


1.584- 


1.595 


Si I 


1.07 


3.70 


0.94 


2.44 


0.98 


2.50 


1.615 - 


1.650 


CO bands 


1.24 


10.5 


0.98 


6.33 


0.95 


9.00 


1.670 - 


1.678 


Al I 


1.32 


3.77 


1.03 


2.14 


1.17 


1.85 


1.70 - 


1.72 


Mgl 


1.56 


4.00 


1.18 


5.57 


0.98 


5.00 



0.854-0.857 


Mg I 


1.02 


1.27 


0.91 


0.51 


0.95 


1.05 


0.865 - 0.869 


Ca I 


1.03 


1.71 


0.89 


1.62 


0.93 


2.15 


1.158 - 1.167 


Fe I 


1.11 


0.08 


0.83 


0.10 


1.09 


0.12 


1.168 - 1.172 


K I 


1.23 


0.04 


0.96 


0.07 


1.11 


0.08 


1.182 - 1.186 


Mgl 


1.33 


0.19 


1.12 


0.09 


1.24 


0.05 


1.187- 1.192 


Fe I 


1.08 


0.09 


0.85 


0.02 


0.90 


0.03 


1.202 - 1.207 


Si I 


1.23 


0.04 


1.05 


0.02 


1.17 


0.03 


1.279 - 1.287 


Ti I 


1.59 


0.23 


1.38 


0.21 


1.59 


0.40 


1.287- 1.292 


Mn I 


0.99 


0.05 


0.88 


0.02 


0.95 


0.03 


1.310 - 1.317 


Al I 


1.23 


0.15 


1.09 


0.10 


1.16 


0.08 
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Table 4. Equivalent Width analysis for Cen X-4 using FS indices 



Feature 


Symbol 


Integration limits 


Equivalent width(A) 


Sil 1.5892 


Wi.59 


1.5870-1.5910 


3.08 ± 0.41 


1^00(6,3) 1.6187 


Wl.62 


1.6175-1.6220 


1.90 ± 0.13 


Nal 2.2076 




2.2053-2.2101 


2.56 ± 0.21 


Fel 2.2263 


Wpel 


2.2248-2.2293 


2.25 ± 0.36 


Fel 2.2387 


WFe2 


2.2367-2.2402 


0.51 ± 0.38 


Cal 2.2636 


Wca 


2.2611-2.2662 


3.30 ± 0.55 


Mgl 2.2814 


WMg 


2.2788-2.2840 


0.89 ± 0.48 


1^00(2,0) 2.2935 


W2.29 


2.2924-2.2977 


4.50 ± 0.44 


i2cO(3,l) 2.3227 


W2.32 


2.3218-2.3272 


3.01 ± 0.43 


^^CO(2,0) 2.3448 


W2.34 


2.3436-2.3491 


1.0 ± 0.59 



No te. — The equiva l ent w idths were calculated using the integration limits 
from iForster Schreibej ||2000|) . The error bar on the measurements are based 
on continuum placement uncertainties. 
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Table 5. Equivalent Width analysis for Cen X-4 using Ali et al. indices. 



Feature 


Symbol 


Integration limits 


Equivalent width{A) 


Nal 2.2076 


WNa 


2.204-2.211 


3.33 ± 0.33 


Cal 2.2636 


Wca 


2.258-2.269 


5.89 ± 1.23 


Mgl 2.2814 


WMg 


2.279-2.285 


1.65 ± 0.57 


^^00(2,0) 2.2935 


W2.29 


2.289-2.302 


5.27 ± 1.09 



Note. — The equivalent widths were calculated using the integration limits 
from Ali, 1995. The error bar on the measurements are based on continuum 
placement uncertainties. 



Table 6. Donor fractions , /, calculated for Cen X-4 with K5 V - MO V field stars 



Band Feature Wavelength range(Atm) / {K5 V) xl f (K7 V) xl f (MO V) xi 



1.97-2.42 


K-band 


1.00 


6.76 


0.90 


7.40 


0.93 


6.76 


1.97-1.99 


Ca I 


1.20 


2.62 


0.85 


3.72 


0.73 


3.83 


2.203 2.214 


Na I 


1.48 


1.39 


1.25 


1.05 


1.15 


1.15 


2.260 - 2.272 


Cal 


1.15 


0.95 


0.98 


1.16 


0.95 


1.17 


2.295 - 2.305 




1.20 


1.38 


1.12 


1.60 


1.00 


1.67 



1.48 - 1.72 


H-band 


0.89 


2.56 


0.93 


2.75 


0.92 


2.76 


1.48 - 1.52 


Mg I 


0.93 


1.24 


0.97 


1.07 


0.98 


1.20 


1.57- 1.582 


Mgl 


0.95 


0.90 


1.04 


1.02 


0.98 


0.85 


1.585 - 1.605 


Si I 


1.10 


1.37 


1.15 


1.60 


1.22 


1.33 


1.615 - 1.65 




0.93 


1.02 


0.93 


0.96 


0.92 


1.20 


1.67-1.68 


Al I 


0.88 


2.46 


0.93 


2.59 


0.90 


2.54 


1.70 - 1.72 


Mgl 


0.71 


2.25 


0.71 


2.21 


0.68 


2.12 





1.182 - 1.186 


Mgl 


0.92 


0.70 


1.03 


0.56 


1.06 


0.72 




1.187-1.192 


Fel 


1.26 


1.01 


1.15 


1.02 


1.23 


1.06 


J 


1.202 - 1.207 


Si I 


1.38 


0.88 


1.71 


1.06 


1.82 


1.28 




1.310 1.317 


Al I 


0.86 


1.55 


0.86 


1.92 


0.80 


0.93 
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Fig. 1. — The NIR spectrum of V404 Cyg. The time-averaged spectrum, shown in gray, was 
obtained after correcting for atmospheric absorption and shifting the individu al exposures to th e 
rest frame of the donor star. The dereddened spectrum, using E(B-V) = 1.303 ()Hynes et al.ll2009l ). 
is shown in black. 
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Fig. 2. — The J-band spectrum of V404 Cyg. Prominent spectral features are labeled at the top of 
the figure and representative error bars are shown near .81^m and 1.35 fiui. 
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Fig. 4. — The K-band spectrum of V404 Cyg with representative error bars. 
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Fig. 5. — The NIR spectrum spectrum of Cen X-4. The time averaged spectrum was obtained by 
correcting the individual exposures for atmospheric absorption and shiftin g them to the res t frame 
of the donor star. The spectrum has been dereddened using E(B-V)=0.1 (jBlair et al.lll984l ). 
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Fig. 6. — The J-band spectrum of Cen X-4. Prominent emission and absorption features are labeled 
at the top of the figure. Representative error bars are shown near wavelengths .94 and 1.25 
/im. 
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Fig. 7. — The H-band spectrum of Cen X-4 with representative error bars. 
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Fig. 8. — The K-band spectrum of Cen X-4 with representative error bars. 
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Fig. 9. — Comparison of V404 Cyg ( shown in black) with a Kl IV field star (shown in gray). The 
template has been normalized to the spectrum of V404 Cyg just blueward of the 2.29//m ^^CO 
absorption bandhead. A zoomed in version of the CO bands (2.29 ~ 2.42 fim) shows that a Kl IV 
star does not match the relative depths of the CO bands in V404 Cyg. 
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Fig. 10. — Comparison of V404 Cyg (shown in black) with a K3 III field star (shown in gray). The 
zoomed in version of the CO bands shows a better match between a K3 III star and V404 Cyg with 
respect to the depths of these absorption features. 
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Fig. 11. — The normalized Kband spectrum of V404 Cyg is shown in black. Overplotted in red is 
a K3 III field star normalized and scaled to / = .97 of K3 III field star. Representative error bars 
are shown in the figure. 
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Fig. 12. — The normalized H-band spectrum of V404 Cyg is shown in black. Overplotted in red is 
a K3 III field star normalized and scaled to / = .98. Representative error bars are shown in the 
figure. 
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Fig. 13. — Comparison of the NIR spectrum of Cen X-4 (shown in black) with a K7 V field star 
(shown in gray). The field star spectrum has been normalized to that of Cen X-4 just blueward of 
i^CO bandhead at 2.29 ^um. 
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Fig. 14. — Comparison of the NIR spectrum of Cen X-4 (shown in black) with a MOV field star 
(shown in gray). The field star spectrum has been normalized to that of Cen X-4 just blueward of 
the ^^CO bandhead at 2.29 /um. 
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Fig. 15. — The normalized H-band spectrum of CenX-4 is shown in black. Overplotted in red is 
a K7V field star, normalized and scaled to / = 0.94. Representative error bars are shown in the 
figure 
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Fig. 16. — The normalized K-band spectrum of CenX-4 is shown in black. Overplotted in red is 
the spectrum of an MO V field star, normalized and scaled to / = 1.00 with representative error 
bars. 
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Fig. 17. — The best f it H-band light curve of V404 Cyg using the orbital ephemeris of 



Casares &: Charles! (|1994| ). T he data points (dia mond) and the equally weighted error bars for 



each point are obtained from ISanwal et al.l (|l996l ) . The solid line shows the best fit to the ellip- 
soidal modulation including 2% dilution from the accretion disk for an inclination of 67° and q = 
16.67 
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Fig. 18. — The best fit H-band light curve of Cen X-4 with zero phase-offset and an inclination of 
35° is shown in solid while the dashed line represents the best fit for an inclination of 36° and a 
phase-offset of .03 . The data points (di amond) and the equally weighted error bars for each point 
are obtained from IShahbaz et al.l (jl993l ). Both include 6% dilution from the accretion disk and are 
plotted for q = 5.00. 



